Genetic variant of the sheep E2F8 gene and its associations with litter size.

The economic efficiency of sheep breeding, aiming to enhance productivity, is a focal point for improvement of sheep breeding. Recent studies highlight the involvement of the Early Region 2 Binding Factor transcription factor 8 (E2F8) gene in female reproduction. Our group's recent genome-wide association study (GWAS) emphasizes the potential impact of the E2F8 gene on prolificacy traits in Australian White sheep (AUW). Herein, the purpose of this study was to assess the correlation of the E2F8 gene with litter size in AUW sheep breed. This work encompassed 659 AUW sheep, subject to genotyping through PCR-based genotyping technology. Furthermore, the results of PCR-based genotyping showed significant associations between the P1-del-32bp bp InDel and the fourth and fifth parities litter size in AUW sheep; the litter size of those with genotype ID were superior compared to those with DD and II genotypes. Thus, these results indicate that the P1-del-32bp InDel within the E2F8 gene can be useful in marker-assisted selection (MAS) in sheep.

Orally Administered Lactobacillus rhamnosus CY12 Alleviates DSS-Induced Colitis in Mice by Restoring the Intestinal Barrier and Inhibiting the TLR4-MyD88-NF-κB Pathway via Intestinal Microbiota Modulation.

Oral ingestion of probiotics is a promising approach to relieving inflammatory disease through regulating the gut microbiota. A newly discovered strain, Lactobacillus rhamnosus CY12 (LCY12), obtained from cattle-yak milk, displayed numerous probiotic properties. These included enhanced viability in low pH and bile environments, adhesion capabilities, and potent antimicrobial effects. The research aimed to explore the beneficial impacts of the novel LCY12 strain on colitis in mice induced by dextran sulfate sodium (DSS) and to elucidate the underlying molecular mechanisms. The results of the study showed that administration of LCY12 effectively helped to reduce the negative effects of DSS-induced body weight loss, disease activity index score, colon length shortening, loss of goblet cells, and overall histopathological scores in the intestines. Simultaneously, LCY12 administration significantly alleviated intestinal inflammation and safeguarded intestinal barrier integrity by enhancing IL-10 levels, while dampening IL-6, IL-1ß, and TNF-α production. Additionally, LCY12 boosted the presence of tight junction proteins. Furthermore, LCY12 hindered the TLR4/MyD88/NF-κB signaling pathway by downregulating TLR4 and MyD88 expression, inactivating phosphorylated IκBα, and preventing translocation of NF-κB p65 from the cytoplasm to the nucleus. The LCY12 also increased specific intestinal microbial communities and short-chain fatty acid (SCFA) production. Altogether, LCY12 oral administration alleviated colitis induced with DSS in mice by improving intestinal barrier function and regulating inflammatory cytokines, SCFA production, and intestinal microbiota.

Copy number variations within fibroblast growth factor 13 gene influence growth traits and alternative splicing in cattle.

Previous researches revealed a copy number variation (CNV) region in the bovine fibroblast growth factor 13 (FGF13) gene. However, its effects remain unknown. This study detected the various copy number types in seven Chinese cattle breeds and analysed their population genetic characteristics and effects on growth traits and transcription levels. Copy number Loss was more frequent in Caoyuan Red cattle and Xianan cattle than in the other breeds. Association analysis between CNV and growth traits of Qinchuan indicated that the CNV was significantly related to chest depth, hip width and hucklebone width (P < 0.05). Additionally, the growth traits of individuals with copy number Loss were significantly inferior to those with copy number Gain or Median (P < 0.05). Besides, we found two splicing isoforms, AS1 and AS2, in FGF13 gene, which resulted from alternative 5' splicing sites of intron 1. These isoforms showed varied expression levels in various tissues. Moreover, CNV was significantly and negatively associated with the mRNA expression of AS1 (r = -0.525, P < 0.05). The CNVs in bovine FGF13 gene negatively regulated growth traits and gene transcription. These observations provide new insights into bovine FGF13 gene, delivering potentially useful information for future Chinese cattle breeding programs.

Exploration of the Polymorphism Distribution of Bovine HMGA2 Gene in Worldwide Breeds and Its Associations with Ovarian Traits.

The high-mobility group AT-hook 2(HMGA2) gene has been widely studied in the context of cancer and animal growth. However, recently, several studies have uncovered its critical role in cell proliferation. A genome-wide association study (GWAS) further suggests that the HMGA2 gene is a candidate gene in fertility, indicating its connection not only to growth traits but also to reproduction, specifically ovarian traits. Thus, this study aimed to analyze the distribution of the HMGA2 gene in 54 bovine breeds worldwide, identify important short fragment variants (indels), and investigate the relationship between HMGA2 and ovarian development. The dataset included genotypic information from a bovine population of 634 individuals (n = 634). After genotyping and analyzing four selected loci, we found that one out of four loci, rs133750033 (P4-D22-bp), was polymorphic. Our results also reveal that this indel of HMGA2 is significantly associated with certain ovarian traits (p < 0.05). Specifically, it has connection with ovarian length (p = 0.004) and ovarian height (p = 0.026) during diestrus. Additionally, we discovered a higher expression of the HMGA2 gene in Asian cattle breeds. In summary, this study suggests that HMGA2 has the potential to serve as an animal fertility testing marker gene. Moreover, these findings contribute to a more promising outlook for the bovine industry.

Transcriptome Analysis of Compensatory Growth and Meat Quality Alteration after Varied Restricted Feeding Conditions in Beef Cattle.

Compensatory growth (CG) is a physiological response that accelerates growth following a period of nutrient limitation, with the potential to improve growth efficiency and meat quality in cattle. However, the underlying molecular mechanisms remain poorly understood. In this study, 60 Huaxi cattle were divided into one ad libitum feeding (ALF) group and two restricted feeding groups (75% restricted, RF75; 50% restricted, RF50) undergoing a short-term restriction period followed by evaluation of CG. Detailed comparisons of growth performance during the experimental period, as well as carcass and meat quality traits, were conducted, complemented by a comprehensive transcriptome analysis of the longissimus dorsi muscle using differential expression analysis, gene set enrichment analysis (GSEA), gene set variation analysis (GSVA), and weighted correlation network analysis (WGCNA). The results showed that irrespective of the restriction degree, the restricted animals exhibited CG, achieving final body weights comparable to the ALF group. Compensating animals showed differences in meat quality traits, such as pH, cooking loss, and fat content, compared to the ALF group. Transcriptomic analysis revealed 57 genes and 31 pathways differentially regulated during CG, covering immune response, acid-lipid metabolism, and protein synthesis. Notably, complement-coagulation-fibrinolytic system synergy was identified as potentially responsible for meat quality optimization in RF75. This study provides novel and valuable genetic insights into the regulatory mechanisms of CG in beef cattle.

A repertoire of single nucleotide polymorphisms (SNPs) of major fecundity BMPR1B gene among 75 sheep breeds worldwide.

The BMPR1B gene is a major determinant of sheep reproductive capacity. Previous studies revealed that Q249R (FecB) is a profound variant of BMPR1B that influences the ovulation rate and litter size in sheep. However, unlike Q249R locus, the full spectrum of single nucleotide polymorphisms (SNPs) within BMPR1B has not been extensively studied. A systematic screen of SNPs in BMPR1B would facilitate the discovery of novel variants that are associated with litter size. This study aimed to investigate SNPs in the BMPR1B gene via whole genome sequence (WGS) data from 2409 individuals of 75 sheep breeds worldwide. Herein, a total of 9688 variants were screened, among which 15 were coding variants and 8 were novel changes. Specifically, we presented the most comprehensive frequency distribution map of the well-known FecB mutation to date. Besides, among the above-mentioned SNPs, one synonymous mutation (g.30050773C > T) was found to be likely under selection and is potentially associated with fecundity in Duolang sheep. Thus, our study greatly expands the variation repertoire of the ovine BMPR1B gene and provides a valuable resource for exploring causative mutations and genetic markers associated with litter size.

Exploring the Sheep MAST4 Gene Variants and Their Associations with Litter Size.

The economic efficiency of sheep breeding can be improved by enhancing sheep productivity. A recent genome-wide association study (GWAS) unveiled the potential impact of the MAST4 gene on prolificacy traits in Australian White sheep (AUW)). Herein, whole-genome sequencing (WGS) data from 26 different sheep breeds worldwide (n = 1507), including diverse meat, wool, milk, or dual-purpose sheep breed types from China, Europe, and Africa, were used. Moreover, polymerase chain reaction (PCR) genotyping of the MAST4 gene polymorphisms in (n = 566) Australian white sheep (AUW) was performed. The 3 identified polymorphisms were not homogeneously distributed across the 26 examined sheep breeds. Findings revealed prevalent polymorphisms (P3-ins-29 bp and P6-del-21 bp) with varying frequencies (0.02 to 0.97) across 26 breeds, while P5-del-24 bp was presented in 24 out of 26 breeds. Interestingly, the frequency of the P3-ins-29 bp variant was markedly higher in Chinese meat or dual-purpose sheep breeds, while the other two variants also showed moderate frequencies in meat breeds. Notably, association analysis indicated that all InDels were associated with AUW sheep litter size (p < 0.05). These results suggest that these InDels within the MAST4 gene could be useful in marker-assisted selection in sheep breeding.

Bovine FRAS1: mRNA Expression Profile, Genetic Variations, and Significant Correlations with Ovarian Morphological Traits, Mature Follicle, and Corpus Luteum.

The amelioration of bovine fertility caused by a multi-factorial problem has always been a hot topic, among which the detection of available target genes is the most crucial. It was hypothesized that the Fraser extracellular matrix complex subunit 1 (FRAS1) gene detected by GWAS is involved in physiological activities such as ovarian development. Herein, unilateral ovaries from 2111 cows were used to examine the mRNA expression profile and polymorphisms of bovine FRAS1 and their associations with fertility-related characteristics. Firstly, it was confirmed that FRAS1 gene transcripts are expressed in various bovine tissues. Then, among five potential insertion-deletion (indel) loci, the 20 bp (named P3-D20-bp) and 15 bp (P4-D15-bp) deletion mutations were confirmed to be polymorphic with linkage equilibrium. Secondly, the P3-D20-bp polymorphism was significantly associated with ovarian weight and corpus luteum diameter in the metaestrus phase and ovarian length in the dioestrum stage. Additionally, both ovarian length and mature follicle diameter in metaestrus are significantly correlated with different genotypes of P4-D15-bp. Thirdly, the transcriptional expression of the FRAS1 gene in groups with a minimum value of ovarian weight or volume was significantly higher than the expression in groups with a maximum value. Instead of that, the more corpus luteum and mature follicles there are, the higher the transcription expression of the FRAS1 gene is. Furthermore, FRAS1 expression in cows with a heterozygous genotype (ID) of P3-D20-bp was significantly higher than others. Eventually, P3-D20-bp deletion could disturb the binding efficiency of WT1-I and Sox2 to FRAS1 sequence according to binding prediction, indicating that mutation may affect gene expression and traits by influencing the binding of transcription factors. Overall, the polymorphisms of P3-D20-bp and P4-D15-bp of the bovine FRAS1 gene significantly correlated to follicle or ovarian traits that could be applied in optimizing female fertility in cow MAS breeding programs.

circUBE3C modulates myoblast development by binding to miR-191 and upregulating the expression of p27.

Noncoding RNAs, including miRNAs (microRNAs) and circRNAs (circular RNA), are crucial regulators of myoblast proliferation and differentiation during muscle development. However, the specific roles and molecular mechanisms of circRNAs in muscle development remain poorly understood. Based on the existing circRNA-miRNA-mRNA network, our study focuses on circUBE3C, exploring its differential expression in fetal and adult muscle tissue of the cattle and investigating its impact on myoblast proliferation, apoptosis, and differentiation. The functional analysis of overexpression plasmids and siRNAs (small interfering RNAs) targeting circUBE3C was comprehensively evaluated by employing an array of advanced assays, encompassing CCK-8 (cell counting kit-8), EdU (5-ethynyl-20-deoxyuridine), flow cytometry, western blot analysis, and RT-qPCR. In vivo investigations indicated that overexpression of circUBE3C impedes the process of skeletal muscle regeneration. Mechanistically, we demonstrated that circUBE3C interacts with miR-191 and alleviates the suppression of p27 through cytoplasmic separation, bioinformatics prediction, dual-luciferase reporter assay, and RIP (RNA immunoprecipitation). Our findings indicate that the novel circRNA circUBE3C competitively binds to miR-191, thereby inhibiting proliferation and promoting apoptosis in bovine primary myoblasts and unveiling a regulatory pathway in bovine skeletal muscle development. These findings expand our understanding of circRNA functions in mammals and provide a basis for further exploration of their role in myogenesis and muscle diseases.

Relationships between the mutations of the goat GATA binding protein 4 gene and growth traits.

Osteogenesis is a complex multilevel process regulated by multiple genes. The GATA binding protein 4 (GATA4) gene has been extensively studied for its pivotal role in bone genesis and bone differentiation. However, its relationship with the growth traits of Shaanbei white cashmere (SBWC) and Guizhou black (GB) goats remains unclear. This work aims to investigate the potential influence of genetic mutations in the GATA4 gene on the growth traits goats. Thus, two Insertion/deletion (InDel) polymorphisms (8-bp-InDel and 9-bp-InDel) were screened and detected in a total of 1161 goats (including 980 SBWC goats and 181 GB goats) using PCR and agarose gel electrophoresis. The analyses revealed that there were two genotypes (ID and DD) for these two loci. In SBWC goats, 8-bp-InDel and 9-bp-InDel loci were significantly associated with heart girth (HG) and hip width (HW). Notably, individuals with DD genotype of 8-bp-InDel locus were superior while those with DD genotype of 9-bp-InDel locus were inferior. Correlation analyses of the four combined genotypes revealed significant associations with cannon circumference (CC), body height (BH), HG and HW. This work provides a foundation for the application of molecular marker-assisted selection (MAS) in goat breeding programs. Furthermore, the findings highlight the potential of the GATA4 gene and its genetic variations as valuable indicators for selecting goats with desirable growth traits.

Global genetic diversity, introgression, and evolutionary adaptation of indicine cattle revealed by whole genome sequencing.

Indicine cattle, also referred to as zebu (Bos taurus indicus), play a central role in pastoral communities across a wide range of agro-ecosystems, from extremely hot semiarid regions to hot humid tropical regions. However, their adaptive genetic changes following their dispersal into East Asia from the Indian subcontinent have remained poorly documented. Here, we characterize their global genetic diversity using high-quality whole-genome sequencing data from 354 indicine cattle of 57 breeds/populations, including major indicine phylogeographic groups worldwide. We reveal their probable migration into East Asia was along a coastal route rather than inland routes and we detected introgression from other bovine species. Genomic regions carrying morphology-, immune-, and heat-tolerance-related genes underwent divergent selection according to Asian agro-ecologies. We identify distinct sets of loci that contain promising candidate variants for adaptation to hot semi-arid and hot humid tropical ecosystems. Our results indicate that the rapid and successful adaptation of East Asian indicine cattle to hot humid environments was promoted by localized introgression from banteng and/or gaur. Our findings provide insights into the history and environmental adaptation of indicine cattle.

The Hoof Color of Australian White Sheep Is Associated with Genetic Variation of the MITF Gene.

Studying the characteristics of mammalian hoof colors is important for genetic improvements in animals. A deeper black hoof color is the standard for breeding purebred Australian White (AUW) sheep and this phenotype could be used as a phenotypic marker of purebred animals. We conducted a genome-wide association study (GWAS) analysis using restriction site associated DNA sequencing (RAD-seq) data from 577 Australian White sheep (black hoof color = 283, grey hoof color = 106, amber hoof color = 186) and performed association analysis utilizing the mixed linear model in EMMAX. The results of GWAS demonstrated that a specific single-nucleotide polymorphism (SNP; g. 33097911G>A) in intron 14 of the microphthalmia-associated transcription factor (MITF) gene was significantly associated with the hoof color in AUW sheep (p = 9.40 × 10-36). The MITF gene plays a key role in the development, differentiation, and functional regulation of melanocytes. Furthermore, the association between this locus and hoof color was validated in a cohort of 212 individuals (black hoof color = 122, grey hoof color = 38, amber hoof color = 52). The results indicated that the hoof color of AUW sheep with GG, AG, and AA genotypes tended to be black, grey, and amber, respectively. This study provided novel insights into hoof color genetics in AUW sheep, enhancing our comprehension of the genetic mechanisms underlying the diverse range of hoof colors. Our results agree with previous studies and provide molecular markers for marker-assisted selection for hoof color in sheep.

Phloretin alleviates palmitic acid-induced oxidative stress in HUVEC cells by suppressing the expression of LncBAG6-AS.

Oxidative stress (OS) is an important trigger of vascular endothelial injury (VEI), which then leads to cardiovascular disease (CVDs). Phloretin was previously investigated to alleviate OS in human umbilical vein endothelial cells (HUVECs) by activating the AMPK/Nrf2 pathway; however, whether phloretin exerts cardiovascular health benefits by targeting non-coding RNAs (ncRNAs) remains unclear. Herein, the whole transcriptome sequencing and lncRNA library building were performed on HUVECs, a commonly used cell line for CVDs study, from different groups in control (CK), palmitic acid (PA, 100 µM), and PA + phloretin (50 µM, G50). KEGG analysis demonstrated that DE-lncRNAs regulated the pathway related to OS and metabolism in HUVECs. LncBAG6-AS was highly expressed under OS stimulation, which was reversed by phloretin co-treatment. Moreover, the MMP, activities of SOD, GSH-Px, T-AOC and GR were significantly ameliorated after interference of LncBAG6-AS, which were consistent with phloretin recover group. Furthermore, the expression of DE-genes from previously reported mRNA sequencing, including MAPK10, PIK3R1, ATP2B4, AKT2, and ADCY9, were significantly changed with LncBAG6-AS interference, indicating that LncBAG6-AS may participate in the process of OS attenuation by phloretin through regulating gene expression. So, the transcriptome sequencing of HUVECs with LncBAG6-AS knockdown was subsequently performed and DE-genes for "NC vs. si-ASO-LncBAG6-AS" were significantly enriched with GO terms, such as apoptosis, response to OS, ferroptosis, and others, which were similar to those observed from KEGG analysis. Overall, this study provides new insights into the molecular mechanisms by which bioactive substances alleviate OS and potential targets for the early prevention and treatment of VEI.

Whole-genome resequencing reveals the genomic diversity and signatures of selection in Romanov sheep.

Romanov sheep are adapted to the extremely cold and harsh environment and display a distinctive grey color. Herein, we analyzed the population structure, genetic diversity, and selection signatures of Romanov sheep based on whole-genome sequencing data of 17 Romanov sheep, 114 individuals from other 10 European breeds. The results of PCA, ADMIXTURE, and NJ-tree showed that the Romanov sheep was closely related to other northern European breeds. A relative high level of genetic diversity, low inbreeding coefficient, and large effective population size was observed in Romanov sheep when compared with other European breeds. We then screened the genomic selection signatures of Romanov sheep using FST, XP-XLP, and XP-EHH methods. The most significant region under selection (CHR14:14.2 to 14.3 Mb) harbored a haplotype that contained MC1R gene. Furthermore, this haplotype was also found in other grey-body breeds including Gotland sheep, Grey Tronder Sheep, and German grey heath sheep, suggesting that it was associated with the unique coat color of these breeds. We also found one region (CHR10:40.8Mb- 41.0Mb) harboring PCDH9 gene which was potentially associated with cold environmental adaptation. In summary, this study identified candidate genes that were associated with the unique grey color and environmental adaptation in Romanov sheep, which provided a basis for understanding the genetic background and utilization of this breed.

Romanov sheep is one of the most famous sheep breeds in the word, characterized by adaptability to harsh environment, high fertility, and unique coat color. Understanding its genetic architecture and signatures is of great value for its conservation and utilization. In this study, we analyzed whole-genome sequences of Romanov sheep as compared with 11 other European breeds, to explore for the population structure, genetic diversity, and selection signatures. We discovered a series of candidate genes that likely play a role in the grey coat color and cold adaptation of the Romanov sheep. In particular, we identified MC1R as a strong candidate gene that determines the grey coat color.

InDel and CNV within the AKAP13 Gene Revealing Strong Associations with Growth Traits in Goat.

A-kinase-anchoring protein 13 (AKAP13) is a member of the AKAP protein family that has been found to be associated with bone formation. Thus, we investigated the AKAP13 gene as a potential candidate gene for molecular-marker-assisted selection (MAS) in breeding. Our aim was to explore genetic variations (InDel and CNV) within the AKAP13 gene of Shaanbei white cashmere (SBWC) goats and analyze their relationship with growth traits. Ultimately, we identified three InDel loci (16-bp deletion, 15-bp insertion, and 25-bp deletion) and three CNVs, and the 16-bp and 15-bp loci were significantly associated with goat body length (p < 0.05). Both the 16-bp deletion variant and the 15-bp insertion variant facilitated an increase in body length in goats. In addition to this, there was a certain superposition effect between 16-bp and 15-bp loci, although there was no linkage. Additionally, the CNV1 locus was significantly correlated with body height and body length of goats (p < 0.05), and CNV2 was significantly correlated with chest depth, chest circumference, and cannon circumference of goats (p < 0.05). Individuals with gain type showed excellent growth performance. In conclusion, the InDel and CNV loci that we have identified could possibly serve as effective molecular markers in goat breeding, which is very essential for improving efficiency and success of breeding. Moreover, our findings provide a new avenue for further research into the function of the AKAP13 gene.

Structural variation and introgression from wild populations in East Asian cattle genomes confer adaptation to local environment.

BACKGROUND: Structural variations (SVs) in individual genomes are major determinants of complex traits, including adaptability to environmental variables. The Mongolian and Hainan cattle breeds in East Asia are of taurine and indicine origins that have evolved to adapt to cold and hot environments, respectively. However, few studies have investigated SVs in East Asian cattle genomes and their roles in environmental adaptation, and little is known about adaptively introgressed SVs in East Asian cattle. RESULTS: In this study, we examine the roles of SVs in the climate adaptation of these two cattle lineages by generating highly contiguous chromosome-scale genome assemblies. Comparison of the two assemblies along with 18 Mongolian and Hainan cattle genomes obtained by long-read sequencing data provides a catalog of 123,898 nonredundant SVs. Several SVs detected from long reads are in exons of genes associated with epidermal differentiation, skin barrier, and bovine tuberculosis resistance. Functional investigations show that a 108-bp exonic insertion in SPN may affect the uptake of Mycobacterium tuberculosis by macrophages, which might contribute to the low susceptibility of Hainan cattle to bovine tuberculosis. Genotyping of 373 whole genomes from 39 breeds identifies 2610 SVs that are differentiated along a "north-south" gradient in China and overlap with 862 related genes that are enriched in pathways related to environmental adaptation. We identify 1457 Chinese indicine-stratified SVs that possibly originate from banteng and are frequent in Chinese indicine cattle. CONCLUSIONS: Our findings highlight the unique contribution of SVs in East Asian cattle to environmental adaptation and disease resistance.

Polymorphism of GTF2A1 gene is associated with litter size in sheep.

General transcription factor IIA subunit 1 (GTF2A1) is required for transcriptional activation from most promoters in eukaryotic transcription. Previous whole-genome association analyses articles have predicted the effect of this gene on lambing in sheep. In the study, a total of nine insertion/deletion (indel) variants named L1 to L9 in the gene were selected to detect in 550 Australian White sheep (AuW) adult ewes. The polymorphisms were detected in four loci (L1, L2, L3 and L8locus) and the polymorphism information content (PIC) values were 0.270, 0.375, 0.372 and 0.314. Further, our study revealed that L1, L2 and L3 loci of GTF2A1 gene were significantly related to the first parity litter size and the polymorphism of L8 was significantly correlated with litter size in the second parity. In details, for the first parity, individuals with the II genotype of the L1 locus had higher little size than that with the ID genotype, individuals with the ID and DD genotype of the L2 locus had higher little size than that with the II genotype, and individuals with the DD genotype of the L3 locus had higher little size than that with the II genotype. All four loci do not follow to Hardy-Weinberg equilibrium, and have no linkage between them. In conclusion, the polymorphisms of GTF2A1 were confirmed and analysis results demonstrated that there are some relevance between difference genotypes and litter size, and these findings may provide new insights for accelerating sheep molecular breeding through molecular marker-assisted selection (MAS).

Goat Pleomorphic Adenoma Gene 1 (PLAG1): mRNA Expression, CNV Detection and Associations with Growth Traits.

The pleomorphic adenoma gene 1 (PLAG1) gene, as the major gene responsible for growth, plays a vital role in myogenesis. Meanwhile, the relationship between copy number variation (CNV) of this gene and growth traits in goats remains unclear. Therefore, this study investigated four aspects: bioinformatics analysis, mRNA expression (n = 6), CNV detection (n = 224), and association analysis. The findings indicated that the gene had a large number of conserved motifs, and the gene expression level was higher in fetal goats than in adult goats. Three CNV loci were selected from the database, among which CNV1 was located in the bidirectional promoter region and was associated with goat growth traits. CNV analysis showed that CNV2 and CNV3 of the PLAG1 gene were associated with growth traits such as body weight, heart girth, height at hip cross, and hip width (p < 0.05), with CNV1 loss genotype being the superior genotype, and CNV2 and CNV3 median and gain genotypes of being superior genotypes. This finding further confirms that the PLAG1 gene is the dominant gene for growth traits, which will serve as theoretical guidance for goat breeding.

A sheep pangenome reveals the spectrum of structural variations and their effects on tail phenotypes.

Structural variations (SVs) are a major contributor to genetic diversity and phenotypic variations, but their prevalence and functions in domestic animals are largely unexplored. Here we generated high-quality genome assemblies for 15 individuals from genetically diverse sheep breeds using Pacific Biosciences (PacBio) high-fidelity sequencing, discovering 130.3 Mb nonreference sequences, from which 588 genes were annotated. A total of 149,158 biallelic insertions/deletions, 6531 divergent alleles, and 14,707 multiallelic variations with precise breakpoints were discovered. The SV spectrum is characterized by an excess of derived insertions compared to deletions (94,422 vs. 33,571), suggesting recent active LINE expansions in sheep. Nearly half of the SVs display low to moderate linkage disequilibrium with surrounding single-nucleotide polymorphisms (SNPs) and most SVs cannot be tagged by SNP probes from the widely used ovine 50K SNP chip. We identified 865 population-stratified SVs including 122 SVs possibly derived in the domestication process among 690 individuals from sheep breeds worldwide. A novel 168-bp insertion in the 5' untranslated region (5' UTR) of HOXB13 is found at high frequency in long-tailed sheep. Further genome-wide association study and gene expression analyses suggest that this mutation is causative for the long-tail trait. In summary, we have developed a panel of high-quality de novo assemblies and present a catalog of structural variations in sheep. Our data capture abundant candidate functional variations that were previously unexplored and provide a fundamental resource for understanding trait biology in sheep.